#Loading spatial data
pacnty.load <- readOGR("./PA_Counties_clip/PA_Counties_clip.shp", layer = "PA_Counties_clip", GDAL1_integer64_policy = TRUE)
## OGR data source with driver: ESRI Shapefile 
## Source: "C:\Users\BadMix\Desktop\CMU Classes\Fall 18\R Shiny\Homework\hw3-atolson\PA_Counties_clip\PA_Counties_clip.shp", layer: "PA_Counties_clip"
## with 67 features
## It has 13 fields
plot(pacnty.load)

#loading economic data
econ.load <- read.csv("./aff_download/ACS_16_5YR_S1903_with_ann.csv")

# merging spatial & econ data
pacnty <- pacnty.load[pacnty.load$NAME %in% econ.load$County,]
pacnty@data <- merge(pacnty@data, econ.load, sort = FALSE, by.x = "NAME", by.y = "County")
#palette defined
pal <- colorNumeric(
  palette = "Blues",
  domain = pacnty$Median_Income)

#leaflet plot with PA median income colored
leaflet(data = pacnty) %>%
  addProviderTiles("Esri.WorldTopoMap", options = providerTileOptions(noWrap = TRUE)) %>%
  addPolygons(color = ~pal(`Median_Income`), popup = ~paste0("<b>", NAME, ":</b> $", format(Median_Income, big.mark = ","))) %>%
  addLegend(position = "topright", pal = pal, values = pacnty$Median_Income, title = "Median Income by County (dollars)")